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SUMMARY 

This study was a theoretical and experimental investigation of thin liquid sheet flows in vacuum. A sheet 
flow created by a narrow slit of width, W, coalesces to a point at a distance, L, as a result of surface tension 
forces acting at the sheet edges. As the flow coalesces the fluid accumulates in the sheet edges. The observed 
triangular shape of the sheet agrees with the calculated triangular res ult. E xperimental results for L/W as a func- 
tion of Weber number, We, agree with the calculated result, LfW = v8 We. The edge cross sectional shape is 
found to oscillate from elliptic to "cigar" like to "peanut" like and then back to elliptic in the flow direction. A 
theoretical one-dimensional model was developed that yielded only elliptic solutions for the edge cross section. 
At the points where the elliptic shapes occur there is agreement between theory and experiment. 
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NOMENCLATURE 


area 

wave propagation velocity 
force 

Froude number 
gravitational acceleration 
length of sheet 
arc length, slit depth 
pressure 

volumetric flow rate 
radius of curvature 

x-position at which edge joins main sheet 
sheet edge shape 


U dimensionless ^-velocity 

u,v,w velocity components in Cartesian coordinate system 


W slit width 

We Weber number 

x ,y,z Cartesian coordinate system 

a shape parameter; ratio of y-z curvature to x-y curvature 

P parameter related to shape parameter 

p viscosity 

£,r|,0 dimensionless Cartesian coordinate system 

p density 

o surface tension 

r sheet thickness 

tp angle between antisymmetric waves and z-axis 

Subscripts 

asy antisymmetric 

c edge cross section 

e edge 

max maximum value 

o initial value 

s surface 

sy symmetric 

x,y,z partial derivatives with respect to, or components in, Cartesian coordinates 

£,r|,0 partial derivatives with respect to, or components in, dimensionless Cartesian coordinates 


INTRODUCTION 

Incompressible, thin sheet flows have been of research interest for many years. Those studies were 
mainly concerned with the stability of the flow in a surrounding gas. Squire (ref. 1) was the first to carry out a 
linear, inviscid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and 
Fraser (ref. 2) did an experimental study of the disintegration of sheet flows using several viscous liquids. They 
also detected the formation of holes in their sheet flows. Hagerty and Shea (ref. 3) carried out an inviscid stabil- 
ity analysis and calculated growth rates for the instability. They compared their calculated growth rates with 
experimental values. Taylor (refs. 4 to 6) studied extensively the stability of thin liquid sheets both theoretically 
and experimentally. He showed that thin sheets in vacuum are stable. Brown (ref. 7) experimentally investigated 
thin liquid sheet flows as a method of application of thin films. Clark and Dombrowski (ref. 8) carried out a 
second order stability analysis for inviscid sheet flows. Lin (ref. 9) introduced viscosity into the linear stability 
analysis of thin sheet flows in vacuum. Mansour and Chigier (ref. 10) conducted an experimental study of the 
break up of a sheet flow surrounded by high-speed air. Lin, Lian, and Creighton (ref. 1 1) did a linear stability 
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analysis that includes viscosity and a surrounding gas. Rangel and Sirignano (ref. 12) carried out both a linear 
and nonlinear, inviscid stability analysis that applies for any density ratio between the sheet liquid and the sur- 
rounding gas. 

Now there is renewed interest in sheet flows because of their possible application as low mass radiating 
surfaces (refs. 13 and 14). Because of their low mass, near immunity to micrometeoroid damage and simplicity, 
sheet flows are excellent candidates for a space radiator system. The objective of this study is to investigate the 
fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern 
the sheet geometry will be derived and compared to experimental results. Since a space radiator will operate in 
a vacuum the analysis does not include any drag force on the sheet flow. 

Thin sheet flows are dominated by the surface tension forces. As a result of surface tension at the edges 
of the sheet, a flow that begins with a dimension, W, perpendicular to the flow direction coalesces to a point at 
a distance, L, in the flow direction. The resulting triangular sheet is ideal for a space radiator. In the following 
analysis results for L/W scaling and variation of sheet thickness t, that include a constant gravity force in the 
flow direction, will be derived. Experimental results for L/W will be compared to the analytical results. Also, 
analytical results for the edge cross-sectional shape will be calculated and compared with experimental results. 


THEORY 

Sheet Thickness and Sheet Length/Sheet Width Scaling 

A sketch of the flow geometry of a thin liquid sheet is shown in figure 1 . Surface tension forces at the 
two edges of the sheet push the edges toward the z-axis. As a result, as the flow moves in the z-direction the 
edge cross-sectional area, A ( „ grows. In order to satisfy mass continuity the edges approach each other and 
finally meet at the point z = L. 

If sheet flow is to be utilized as a space radiator then the dependence of the sheet length to initial width 
ratio, L/W, and sheet thickness, r, on the flow conditions must be known. In this section the L/W dependence on 
the Weber number, We, will be determined. Also, since all experiments have been performed in the Earth’s 
gravity field the L/W and t dependence on a constant gravity force will be included. 

First consider the total force exerted by the surface tension on the edges of the sheet. Consider a cross 
section of the sheet taken perpendicular to the flow direction as shown in figures 1 and 2. The surface tension 
pressure is o/R c , where a is the surface tension of the fluid and R c is the radius of curvature (ref. 15). Referring 
to figure 2, the surface tension force on the infinitesimal surface area dA = dl dz is 

dF=- — dldz (1) 

R, 


The component of this force in the x— direction is 


dF = - — sin 4 ) dl dz = ~ — ds dz 
X Rc Rc 


( 2 ) 
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where 5 is the cross-section shape and the radius of curvature, R c , is given by 


dh 

1 dx 2 

X 


In using (3), any curvature in the y-z plane is being neglected compared to the curvature in the x-y plane. This 
is a good assumption for long {LfW > 1) sheet flows. 

Substituting (3) in (2) and then integrating over the surface from s = 0 to s = t/2 will yield the x-direction 
force on the upper half of the sheet edge. Therefore, the total edge force in the x-direction per unit length in the 
z-direction is the following: 


1 + 

r ds V 


X) . 



(4) 


where the prime denotes differentiation with respect to x. Making the substitution u = s' in (4) yields 



du 


(5) 


so that 




dz 


= 2a 


( 6 ) 


The significance of this result is that it is independent of the cross-sectional shape of the sheet. It applies as 
long as the edge cross section has infinite slope (s" — *• °°) at the edge (x = 0) and zero slope (s' = 0) where it 
joins the constant thickness region (x = r). 

The x-direction force given by (6) can be used in the integrated momentum equation for the control 
volume shown in figure 2. This coordinate system is moving with the edge velocity, u g . Therefore, the x-direction 
velocity is u = -u e at x = r and u = 0 at x = 0. The integrated momentum equation combined with (6) yields the 
following result. 


2 dF x 

pt u = — — = 2 a 
dz 


(7) 
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Therefore, 


u 


e 


2a 


N PT 


( 8 ) 


The velocity, u g , is the velocity that a free edge of a thin liquid sheet will move regardless of the edge 
shape. Taylor (ref. 6) derived this result earlier without considering the cross-sectional shape of the edge. This is 
also the phase velocity, y , of the antisymmetrical wave that can exist on a plane sheet (refs. 5 and 6). Thus 
for a constant thickness sheet moving with velocity w in the z-direction (fig. 1) stationary antisymmetrical waves 
will make an angle, tp', with the flow direction defined by the following relation (ref. 5) 


sin ip' = 



w w 


( 9 ) 


Also, referring to figure 1, 


u e 

tan ip = — 


w 


( 10 ) 


Therefore, sin ip' = tan ip. For long sheets (LfW >1), tan ip » sin ip, so the sheet edge and stationary antisym- 
metrical waves (lines of constant phase) will be approximately parallel. Experimentally, these waves are 
observed (fig. 16) as they have also been observed in other experiments (refs. 3, 6, and 10). 

Equation (8) gives the sheet edge velocity, u g . With this result and the continuity equation, analytical 
expressions for t(z), x(z), and UW can now be developed. Referring to figure 1(a), at z = 0 the flow will be 
moving in the z-direction with initial velocity, w () . Eventually, fluid that begins at some position, x, at z = 0 will 
be captured by the edge with cross-sectional area A g . Therefore, for steady-state conditions the following conti- 
nuity equation must be satisfied. 


A c w g + T xw = 


1 

2 




(ID 


In the edge cylinder the z-direction velocity is w g , while in the main sheet the z-direction velocity is w. The w 
velocity will be determined by the gravitational force, whereas w g is determined by the combined gravitational 
and surface tension forces. Rather than use the complete continuity equation (11), the following simplifying 
approximation is made. We assume the flow velocity remains entirely in the z-direction until it reaches the 
edge. Therefore, 

T XW = X 0 XW 0 (12) 

Near the slit that creates the sheet flow viscous forces will decelerate the flow. However, Brown (ref. 7) 
has shown this effect is only important near the slit. Thus, if the gravity force is in the z-direction (fig. 1) the 
main sheet flow velocity, w, is governed by the following equation of motion. 
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dw dw 

— = w = g 

dt dz 


( 13 ) 


where g is the gravitational acceleration. Integration of (13) yields the following. 

w 2 = w 1 + 2 gz 


(14) 


Also, using the sheet edge velocity in the ^-direction given by (8), 


dx dx 

— = w — = -u. = - 
dt dz 


N 


2o 

pr 


(15) 


Now substitute (14) and (12) into (15) to obtain 


dx 

dz 


= -j2Wew c 


w n + 2 gz 


- 1/4 


(16) 


where the Weber number, We, which is the ratio of the dynamic pressure to the surface tension pressure, is 
given by 




(17) 


Integration of (16) yields the following result for the sheet shape. 


— = 1 - - Fr \[mre 
W 3 


1 + — — 
Fr W 


\ 3/4 


- 1 


(18) 


where the Froude number Fr, which is the ratio of the kinetic energy to the gravitational energy, is given by 

2 

Fr = — 


(19) 


When x = 0 the sheet flow will have reached its maximum length, z = L. Therefore, setting x = 0 in (18) yields 
the following result for L/W. 


_L 

W 


El 

~2 


MI3 


i + !JL 


L V 


2 Fr N 


We 

T 


( 20 ) 


- 1 
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The importance of We in determining UW becomes obvious when (20) is expanded for large Fr. For all the 
experiments carried out to date Fr > 10. For the case of no gravity field (g = 0), Fr and (20) yields 



We 


\ 8 


( 21 ) 


Since JWe ~ w Q , (21) indicates that L will be a linear function of w 0 for fixed IV and t q . This was observed in 
the experiments of Mansour and Chigier (ref. 10) and Chubb and White (ref. 13). Equation (21) for L/W and the 
complete expression (20) for several values of Fr are shown in figure 3. For Fr > 20 there is only a small devi- 
ation from the Fr -*■ 00 (g = 0) result. Thus for the experimental conditions of this study, (21) can be used to 
estimate L/W. 

Using the result (21) for L/W in (18) yields the following equations for the sheet shape in terms of the 
dimensionless coordinates 2 x/W and zJL. 


where 


2x 

~W 


= 1 - a 


1 + 


>-V 



( 22 ) 


a s iFr^S/We 


(23) 


and 


( 

p - 1 

V 



- 1 


(24) 


The quantity, a, can be called the sheet shape parameter since it determines the sheet shape. For most cases of 
experimental interest a > 1. For a — (22) becomes 


— = 1 - — for a -*• 00 (25) 

W L 


Therefore, the sheet is triangular in shape for a » 1. The sheet shape as given by (22) is shown in figure 4(a) 
for several values of the shape parameter. For a > 1 the sheet is essentially triangular in shape. Indeed, this is 
the result observed experimentally. 

Now consider the variation of sheet thickness, r(z). Combining equations (12), (14), and (20) yields the 
following result. 


1 + P- 


V 1/2 


(26) 
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Similar to x(z), the sheet thickness variation t(z) depends only on the shape parameter, a. Also, for no gravity 
force, (3 = 0, and (26) yields r/0 o = 1; in this case the sheet thickness is a constant. The sheet thickness calcu- 
lated using (26) for several values of a is shown in figure 4(b). These theoretical results show that significant 
thinning of the sheet occurs even for a > 1. However, for a > 10 negligible thinning occurs. These results have 
not been verified experimentally. 


Cross-Sectional Shape of Thin Liquid Sheet 

In the preceding section, results for L/W and r were derived without the need for knowing the cross- 
sectional shape of the sheet edge. Now the analysis for the cross-sectional shape ( s(x £ ) in fig. 2) will be pre- 
sented. The cross-sectional shape of the sheet edge is important in evaluating the stability of the sheet. The large 
curvature (small radius of curvature) that exists in the region where the edge joins the sheet (x = r in fig. 2) can 
lead to instability. This problem is currently being investigated. 

As pointed out in the last section, the effect of the z-direction gravity force is negligible for the sheet 
flows of interest. Therefore, only the surface tension force must be considered. In equation (3) of the previous 
section, the curvature of the sheet in the x-z plane was neglected so that the surface tension force was confined 
to the x-y plane. This is a reasonable assumption for flows of sufficient length ( L/W > 1). In that case, the 
cross-sectional shape, s(x,z), increases slowly in the z-direction so that » s u and only curvature in the x-y 
plane is important. The radius of curvature is then given by (3). 

For the most general case, however, curvature in the x-z plane must be included. In the following 
analysis, the velocity field is assumed to be u = [u{x,z),v{x,yj),w ( J{, that is, the flow velocity is constant in the 
z-direction and the x-component of velocity is independent of y. Assuming irrotational flow, Bernoulli’s equa- 
tion applied at the surface yields: 



_L + _L 






a 



(27) 


where the subscript s denotes conditions at the sheet-vacuum boundary. Since the flows of interest are in a 
vacuum the only pressure appearing in (27) is the surface tension pressure, afl//?^ + 1//^], where 


l_ + ± 


R 


xy 
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d 2 s + d 2 s 
dx 2 dz 2 


T3/2 



(28) 


The right-hand side of (27) is the surface tension pressure that occurs at x = 0 in figure 2 and is a function of z. 


1 



*yz 


’*_lx=0 


(29) 
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The sheet surface boundary condition for steady state conditions (A-4) is the following: 


ds L ds 
v — u + w 

5 s dx ° dz 


( 30 ) 


Substituting (30) in (27) we obtain 


PT 

2a 


“A 1 + s x) + 2u s w <?x s z + w o s z 


- T. 


s xc + 


2 2 \ 3/2 ^ 

1 + ** + s zl 


(31) 


where the subscripts x and z denote partial differentiation. The edge velocity, u g , is given by (8) and the Weber 
number, We, by (17); substituting these in (31) gives: 


^ * .J) * ^ 

u e 


u g ^ 2 


We We 2 

S r S. + s , 

X z 2 z 


- T - 


+ ^ + T 


2 2\ 3/2 R e 

l+s x + s z) 


(32) 


If the following dimensionless variables are defined, 


lr - _ = ^ n _ Z — rr s 

(J ■ 5 tj ) v ■ ^ r 9 » - - 

r /2 t /2 t /2 t /2 5 


(33) 


then (32) becomes: 


+ ^60 “ -2 


t/: 



( We \ 2 t 

i 2 2 

+ 

W ne " 

1 + q ? + q e 


3/2 


(34) 


An expression for velocity, I/ s , is obtained from the continuity equation derived in Appendix A. By the 
assumption u = m(x,z), we have u s = u, where u is the average velocity in the x-direction. Equation (A-7) yields 
the following: 


U. 


1 dA/dz 
q dA c /dz 


where 


A(x,z) = 2 f x s dx' 
•JO 


(35) 


(A-8) 
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and 


A c (z) = 2 J^ z) s dx (A-9) 

is the total cross-sectional area of the sheet edge. Equation (34) is the equation for the sheet edge shape, T|, as a 
function of the dimensionless coordinates Z, and 0. To simplify the solution of (34), the 0 dependence of the 
cross-sectional area. A, is assumed to occur through the edge length dimension, r{z ). In other words, in rewriting 
(A-8), 



it is assumed that the integral in (36) is independent of z- In other words, the sheet edge cross-sectional shapes 
at each z location are similar. Therefore, 


where 


and 




1 dAJdz 
t) dAJdz 


1 5 A/50 _ 1 A 

tf dAjdQ "n T c 



(37) 


(38a) 


(38b) 


Since A depends only on Z, equation (37) shows that U s is now dependent on 0 (or z) only through A c . 
Substituting (37) into (34) gives the following result for r|(l;,0): 


_ \2 


955 + ^ee - 


Jc 

J 




T 

~Rl 


2 2 

1 + + r]Q 


3/2 


(39) 


It is desirable to eliminate the 0 dependence in (39) so that an ordinary differential equation can be obtained. 
This is done by making the following assumptions. First, it is assumed that in the region of the sheet edge, 
r|c » t)q. This is a good assumption for long sheets (L/W > 1) where the edge area, A c , grows slowly. Then 
(39) becomes: 
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3/2 


(40) 


1 

l 

( \ 

A 

2 L 2\ T 

1 + q 1 

^ ’lee - 2 

h 2 

w 

' T e \ 


In solving (40), an approximation will be made to treat q 00 as an independent parameter. 

As can be seen, by neglecting the q 0 terms the only dependence on the fluid properties (p,o), occurs 
through A c , which depends on W/L (eq. (21)). Thus, for a given value of A c the sheet cross section is the same 
for all fluids. Solution of equation (40) depends on the parameters A c , T/R e , and q^, plus the boundary condi- 
tions at £ = 0 and 2; = r. 


At £ = 0: q = 0, —1 -► ~ (42a) 

At £ = nq = l, 43- = ° (42b) 

The total dimensionless cross-sectional area, A c , must satisfy the condition A c = A (r), so that it is not a free 
parameter. Instead A c can be calculated directly from the flow conditions using the continuity equation (11) 
with the Fr —< ► 00 assumption so that w = w Q = w e and t = t q . Therefore, combining equations (11) and (25) 
yields: 


A 


c 


W t 


o Z 


2 L 


or 


A 


c 


W z 

~o T 


(42a) 


(43b) 


As already mentioned, q 00 is small for the flows of interest ( UW > 1). However, the two sets of bound- 
ary conditions imply that two independent parameters, namely t IR e and q^, are needed in order to obtain a 
physically meaningful solution. If q 00 is neglected in (40), it is found numerically that a physically meaningful 
solution can be obtained for only a unique value of A c Since q 00 cannot be neglected, a simplifying assump- 
tion is required to reduce (40) to an ordinary differential equation. The simplest assumption is that q 00 is a con- 
stant. The procedure for solving (40) using this assumption will now be described, after which it will be shown 
that this assumption does not yield the true solution. A different assumption will be made, namely = aqcc. 
Solutions were investigated using the fourth order Runge-Kutta numerical method (ref. 16). 

At Z, = 0, q = 0, and q^ is infinite, so (40) must be solved analytically near this point. For £ — * 0, A ~ 
q£ (eq. (38a)), so that A /q -*■ 0 and (40) can be approximated by 



-T 

! + q^ 

3/2 2 R e 


for £ ~ 0 


(43) 
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This is the equation for a circle with radius 2R g It and centered at (£,q) = (RJrfi). Therefore, the following are 
used to calculate t) and q^ at the first step in the Runge-Kutta integration of (40): 




r 



AR 

q - 


e - 5 

N 


T 


q ? * 


2R e /r - ? 

0 


(44) 


(45) 


The iterative method for solving (40) consists of setting initial values for T IR e , q 00 , A c , and r, then run- 
ning the Runge-Kutta routine to get revised values for A c and r . Eventually, A c and r converge to a set of val- 
ues. Since A c is kept equal to A (r), and r is defined where = 0, the value of q must be checked to see if it 
equals 1 at r. If not, a new guess for q ee is made. The process is repeated until a physically meaningful solution 
is obtained for the value of T IR e assumed. 

This method was applied for several values of T IR e , and it was determined that a unique solution exists 
for all values of A c . The largest value of A c investigated was 325 200, and the largest A c which ^arises in the 
experiments conducted to date is 2350. Thus it is concluded that a solution to (42) exists for all A c of interest, 
and hence for all sheet flow conditions. Figure 5 shows a typical sheet edge shape for A c = 1 176, assuming q 00 
is a constant. 

Under the assumption that q^ is a constant, equation (40) breaks down for A c < 457. Figure 6 shows the 
two characteristic dimensions of the edge shape solution r and q,^ (fig. 5) as a function of A c (=[z/t][WZL]). 
Both quantities appear to follow a trend for A > 457, but that trend breaks down for lower areas. Note that 
A = 457 is the area at which q 00 = 0. 

An examination of the behavior of qgg as a function of A c suggests that a constant q 00 is not a good 
assumption. First, it is known that the curvature in the thin part of the sheet must be zero, since this is essen- 
tially a planar region of the sheet. But at the sheet edge, we have assumed q 00 constant. Except for the unique 
position on the sheet where q 00 = 0, there is a discrepancy. Second, by figure 6, we see that q 00 < 0 for small 
edge areas, but positive curvature would be expected intuitively. The fluid which is just leaving the slit (where 
z/L = 0 and hence A = 0) is planar, so q 0 = 0. As the fluid flows down the z-axis, the edges grow, so we 
know q 0 > 0. Therefore, we must have q 00 > 0 for at least the lowest range of A c . Since, from figure 6, the 
edge shape solutions break down where the z-curvature is negative, this intuitive violation appears to be related 
to the breakdown in the theory. Finally, we have assumed the edge shapes to be similar at all stations along the 
sheet in order to obtain equation (37). This implies (q,^/^) would be approximately constant, but it is clear 
from figure 6 that such is not the case. The assumption that q 00 is constant therefore produces an unsatisfactory 
set of solutions. 

The remedy is to make a better assumption on q 00 which avoids these violations. Since surface tension is 
presumed to act uniformly in all directions, a suitable assumption is to consider the curvature in the y-z plane to 
be related to the curvature in the x-y plane. This is done by setting q^ = a q^, where a is a constant which 
replaces q 00 as a free parameter. Doing this, (40) becomes 


% = 




r _ \ 2 , v 

l 

1 

A 

1 + qj - J_ 

2(1 + a) 

h 2 

A, 

1 W 



(46) 


12 


Near ^ » 0, this becomes 




2 

3/2 

1 + q 5 



-T 

2R e (\ + a) 


for 5 “ 0 


(47) 


which is the equation for a circle with radius 2(1 + (X)RJt and centered at the point (£,q) = [(1 + a)R e /t,0]. 
Hence for ^ « 0 we have: 


tl 




_ 



4(1 + a)R e 


5 

- K 



T 




2(1 + a)R e /x - 5 
tl 


(48) 


(49) 


Again, the Runge-Kutta algorithm reveals that solutions exist for all A c Figure 7 shows a typical cross-sectional 
shape for an A c of 1231. Note that the edge shape is changed considerably as compared with the shape of fig- 
ure 5. In figure 8, r, q maY , rq maY , RJt and a are plotted as a function of A c In figure 8(b) the product rtlmax 
is shown as a function of A c . This linear dependence is the following. 

= 0-731 rq^ (50) 


If the shape is purely elliptical then using equation (38b), 


A 


c 





max 


0.785 rq 


max 


(51) 


Comparing equations (50) and (51) we see that the cross-sectional shape is nearly elliptical. 

As figure 8 shows, the assumption that q 0 Q = ccq^ does not result in any of the violations which occur- 
red when r|QQ was assumed constant. First, r and q max are in approximately the same ratio for all A c , indicating 
similarity of the sheet edge shape along the flow. This is consistent with the simplifying assumption made about 
the z dependence of A, (eq. (37)). Furthermore, the quantity (filmax) is directly proportional to A c , as expected. 
Since q^ = 0 on the thin portion of the sheet, we have q 6 Q = 0 in this region the expected result. Finally, over 
most of the edge we have q^ < 0. Since for small A c , a < 0, we have q Q Q > 0 for small edge areas. By the 
preceding discussion, this is the expected result. Note that at A c = 0, we must have r = 0, and therefore q is 
discontinuous at ^ = 0. Thus the theory cannot predict the shape solution as A c -* 0. 

The possible existence of oscillatory solutions (i.e., solutions having more than one local minimum for 
5 > 0) was investigated numerically, and it was determined that such solutions do not exist. Instead, 
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equation (40) quickly diverges for % > r. Occasionally oscillatory solutions were observed, but these occurred 
only when the assumed value of A c was much larger than its true value obtained by integrating the area under 
the shape curve. Also, such solutions were found to diverge when the step size was reduced. Therefore, only 
one solution, such as that shown in figure 7, appears to exist for each A c . 


EXPERIMENTAL RESULTS 

The objectives of the experimental program are twofold. The first is to verify the scaling relation (L/W 
versus We from (21)). The second is to measure the edge cross section ( r and T] max ). In the previous studies of 
Chubb and White (ref. 13), and Chubb and Calfo (ref. 14), a limited amount of scaling data was obtained. Most 
of those results were obtained with small slit widths (W < 3.5 cm). For a space radiator system, however, larger 
slit widths will be used. Therefore, the results reported here are for larger slit widths (12.5 cm < W < 25 cm). 
All the L/W versus We data were taken using low vapor pressure Dow Coming 705 silicone oil at vacuum con- 
ditions (less than 1(T 2 torr). However, the edge cross-section results were done in the atmosphere using water as 
the working fluid. Also, the edge shape data were obtained using small slit widths. 


Experimental L/W Scaling Results 

A schematic of the experimental facility used for the L/W scaling results is shown in figure 9. It consists 
of a 30-cm inner diameter stainless steel pipe 3.5 m long. The axis of the pipe is aligned with the gravity field. 
Vacuum conditions exist inside the pipe with the pressure less than 1(T 2 torr. Flow of the Dow Coming 705 sil- 
icone oil through the narrow slits is maintained by pressurizing a 40-gal reservoir with nitrogen (nitrogen gas is 
separated from the oil by a diaphragm to prevent gas entering the oil). The design of the slits is shown in fig- 
ure 10. A length to slit thickness ratio, lh 0 = 7 was chosen to assure that fully developed laminar flow would 
result. 


Four slit widths were used: W = 12.5, 18.75, 23.5, and 25 cm. For W = 12.5 cm, three thicknesses, 

T 0 = 100, 200, and 300 pm, were run. For W = 18.75 cm thicknesses of t q = 200 and 300 pm were run. Thick- 
nesses of t 0 = 100 and 200 pm were run for W = 23.5 cm, and for the W = 25 cm slit a single thickness of T 0 
= 300 pm was investigated. The volumetric flow rate, Q, through these slits was found to satisfy the following 
relation (ref. 7). 


AP tJw (52) 

12p/ 


Where A p is the pressure drop across the slit and p is the viscosity coefficient. This result was derived for 
steady flow of liquid between parallel planes by Lamb (ref. 17). 

The flow velocity, w Q , was determined from the measured volumetric flow rate, Q. 


w 


o 



(53) 


Flow rates were determined by measuring the pressure drop across a calibrated orifice in the supply line to the 
slits rather than using (52). The velocity, w Q , was used to calculate the Weber number, We, using (17). The 
length, L, of the sheet was measured with the aid of a video camera mounted on a translating table shown in 


14 


figure 9. It is not possible to photograph a complete sheet in this facility. However, figure 1 1 shows a complete 
sheet issuing from a small slit (W = 3.4 cm). This photograph is from Chubb and White (ref. 13). Note the tri- 
angular shape of the sheet, which agrees with the calculated result (fig. 4(a)). 


Scaling results are shown in figure 12 where L/W is plotted as a function of We. As can be seen the 
experimental data points all lie close to the calculated result (21). The properties of Dow Coming 705 oil are 
temperature-dependent with the density and surface tension given by the following relations that were obtained 
from the Dow-Coming Company of Midland, Michigan. 


p = 169 - 1.625T fpl = kg/m 3 , [T] = K 


(54a) 


a = 0.08614 - — _ — 
6502 


[d] = N/m, [T] = K 


(54b) 


For the experimental results shown in figure 12, T = 98 °C. 

The results in figure 12 are for vacuum conditions. However, experiments in air using water as the sheet 
fluid yield results that agree with (21). Only a limited range of Weber numbers (50 < We < 300) was investi- 
gated in the atmospheric experiments. For high velocity flows in air an instability will occur that breaks up the 
sheet flow. 


Cross-Sectional Shape Experimental Results 

The sheet was formed in the atmosphere using water as the working fluid. The photographic technique 
employed is shown in figure 13. The dashed line represents the line of sight of the camera. A beam splitter and 
two front faced mirrors were used in order to view both the front and side views of the sheet on the same 
photograph. The tip of a wire was placed near the sheet edge in order to insure that both images on the photo- 
graph were at the same point on the cylinder. Because the sheet edge moves away from the camera in the side 
view, there was only one place in this view that was in focus. As a result, the camera was set on the tripod with 
the ability to move forward and back. If the focus on the camera is not adjusted, then any place that is in focus 
will have the same degree of magnification as any other place that is in focus. Realizing this, it was possible to 
take a picture of the front view in focus, tape the lens to insure a constant focal length, and move the camera 
forwards and back in order to get a number of points along the edge view of the sheet in focus. 

The camera used was a Nikon F4 equipped with a fully extended bellows and a 200-mm Nikon lens, 
with a +2 diopter close-up lens. Kodak TMAX 400 black and white film was used. Three different slit thick- 
nesses (50, 80, and 100 pm) were used, and one set of data was collected from each. All the slits were 3.4 cm 
wide. 


Once a full, steady liquid sheet was formed the length of the sheet was measured. The wire tip was 
placed near the edge, close to the top of the sheet and the distance from the tip of the wire to the top of the 
sheet was measured. The camera and mirrors were placed so the wire would have the same length in each view 
near the bottom of the frame. Two pictures were taken showing both views, and the lens was fixed in place to 
ensure an equal degree of magnification for each picture. Two pictures were then taken showing just the front 
view for greater clarity. Four or five pictures were taken showing different points in the edge view in focus. 
This procedure yielded four or five data points along approximately 1 in. of the edge. 
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At this point, a picture of a millimeter scale was taken in order to determine the degree of magnification 
of the pictures after developing. Then the wire was moved down allowing the top of the next set of photographs 
to be at the same point as the bottom of the previous set. The height of the wire tip was measured again and the 
next set of pictures taken. A picture of the millimeter scale was taken after each set of pictures to check that the 
degree of magnification was not changing. The length of the sheet was periodically measured to check that the 
length of the sheet was staying constant. 

Enlarged photographs were made for taking the r and measurements. The pictures of the scale were 
measured with a millimeter scale to determine the degree of magnification. Knowing the position of the wire 
tip, the blown-up photograph of the scale could be used to determine the height of the point that appeared to be 
in focus in each of the edge shots. The cylinder widths and thicknesses were measured to a resolution of a 64th 
of an inch. 

Knowing the dimensions of the slit, the total length of the sheet, and the height at which the data points 
were taken, the dimensionless cross-sectional area of the edge, A c , could be calculated using (42b). Knowing 
the degree of magnification and the dimensions of the slit, q and r could be calculated using (33). These 
data and the theoretical predictions are shown in figure 14. 

Figure 14 reveals a very interesting phenomenon. The edge appears to be oscillating in the flow direc- 
tion, z. The edge quickly flattens out and then reforms in its basically elliptic shape. When the edge is in its 
elliptic shape it is in good agreement with the theoretical predictions. The details of this cycle are shown in 
figures 15 and 16. Figure 15 shows a cross section of the edge in the x-y plane; figure 16 is a photograph taken, 
as described in figdure 13, of edge and plan views of the sheet. Figure 16(a) is a photograph of the edge very 
near the top of the sheet, where the edge takes on a basically elliptic shape. Figure 16(b) shows the section of 
the edge just below that shown in figure 16(a), where the edge goes through the flattening out and reforming 
process. At the bottom of each photograph the front and side views of the wire tip can be seen. Also notice the 
antisymmetrical waves (with lines of constant phase running parallel to the sheet edge) that appear next to the 
edge in the plan view in figure 16. 

The edge begins in its nearly elliptic shape, shown in figure 16(a). Due to the high curvature where the 
edge connects to the sheet, the edge quickly flattens out to a "cigar" like shape, shown near the top of figure 16(b). 
This appears in figure 14 as an increase in r and a decrease in r] max . Instead of simply regaining its elliptic 
shape, the edge goes through a rebuilding process where the edge has a "peanut" like shape. As the area of the 
edge increases, q max increases first on the inside of the previous edge. This can be seen best a third of the way 
down the edge view of figure 16(b), where the inner part of the edge can be seen behind the outer part of the 
edge. It is not clear whether the edge necks down between these two parts. This inner part of the edge grows 
and engulfs the outer part to return to its elliptic shape. This appears in figure 14 as an increase in q max and a 
decrease in r. This is most noticeable in the 100-|im slit. 

As discussed earlier, the only edge shape solutions obtained analytically are elliptic (fig. 7). No "peanut" 
or "cigar" like shapes, which were observed experimentally (fig. 15), were calculated. Assumptions were made 
about the z-direction derivatives, t|q and r|gg. in order to reduce the two-dimensional problem (34) to a one- 
dimensional problem (46). Also, in (36) the edge shapes were assumed similar at each z-location in order to 
reduce the problem to one dimension. However, the "peanut" and "cigar" shapes observed indicate that the 
assumptions q 0 « q^, t|gg = aq^, and similarity of the edge cylinders do not apply in those cases. Solutions 
like the experimental "cigar" and "peanut" edge shapes require a three-dimensional formulation. 
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SUMMARY 


This study was a theoretical and experimental investigation of the fluid mechanics of thin liquid sheet 
flows in vacuum. Surface tension forces at the edges of the sheet result in a triangular shape for the sheet of 
width, W, and height, L. The flow coalesces to a point at the distance, L, with the fluid being accumulated at the 
sheet edges. Theoretical results that include the gravity force acting in the flow direction yield a triangular shape 
for the sheet. Extensive experimental d ata fo r LfW as a function of Weber number. We, agree with the theoreti- 
cal result for no gravity force, LfW = / We/S. Theoretical results show that for Froude numbers of experimental 
interest the gravity force has negligible effect on LJW. 

Solutions for the edge cross-sectional shape were obtained from a one-dimensional model that includes 
only the surface tension forces. All calculated edge cylinder solutions are elliptic in shape. However, experimen- 
tally it was found that the edges oscillate in the flow direction between elliptic, "cigar" and "peanut" like cross- 
sectional shapes. There is agreement between the theory and experiment at the points where the elliptic shapes 
occur. To account for "cigar" and "peanut" like shapes requires the theoretical model to be extended to three 
dimensions. 
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APPENDIX A 


Continuity Equation Solution 

The continuity equation for steady-state incompressible flow is the following. 

du + dv + dw _ q 

dx dy dz 


(A-l) 


where u, v, and w are respectively the x-, y-, and z-direction velocities. Referring to figure 2, if (A-l) is inte- 
grated from 0 < y < s(xg) then (A-l) becomes the following. 


— ( su) + —( sw) + v - u 
dx J dz 5 


s 


ds 

dx 


ds f, 
w s — = 0 

s dz 


(A-2) 


where u and vv are average velocities defined as follows: 



(A-3) 


and since the flow is symmetric about the *-axis, v(y = 0) = 0. 

At the sheet surface the fluid y-direction velocity, v 5 , must equal the surface speed, ds/dt. Therefore, for 
steady-state conditions. 


ds _ ds 
dt s dx 


+ 


ds 

s dz 


(A-4) 


Using (A-4) in (A-2) yields the following 

dM + d(sw) _ Q (A-5) 

dx dz 

Since the gravity force and surface tension force (curvature in the y-z plane) in the z-direction are being 
neglected, w z « m x . Therefore, 

d(su) = _-ds_ (A-6) 

dx dz 

Beginning at x = 0 where s = 0, and ii = 0, integrate with respect to x to obtain: 

su = - _ vv (A-7) 

2 dz 
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where 


A(x,z) = 2 j^ x s{x! £)dx! (A-8) 

is the sheet edge cross-sectional area for 0 < ^ <x (fig. 2). When jc = r(z), we have A = A c and sit = -(t/2 )u e , 
where A c is the total sheet edge cross-sectional area. Therefore, 

A c (z) = 2 J^ z) s(x' r z)dx' (A-9) 
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(a) Top view. 






Figure 2. — Sheet edge control volume in coordinate system moving 
with edge velocity, Uq. 



Figure 4. — Sheet shape as a function of sheet shape parameter, a. 
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Figure 5. — Sample calculated cross section. Assumes 
2 - curvature constant. 



(a) 



<b) 

Figure 6. — Dependence of sheet edge parameters on edge 
cross-sectional area. Assumes y-z-curvature constant. 



Figure 7. — Sample calculated cross section. Assumes y-z 
curvature proportional to x-y curvature. 
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(c) Product of dimensionless edge dimensions. (d) Ratio of y-z plane curvature to x-y plane curvature. 

Figure 8. — Dependence of sheet edge parameters on edge cross-sectional area. Assumes y-z curvature propotional to x-y curvature. 
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Figure 9. — Schematic of experimental facility for determining 
UW scaling. 
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30 




Figure 10. — Typical slit plate for producing sheet flow. Figure 1 2. — Comparison of theoretical and experimental sheet 

scaling. 




Figure 13. — Photographic set up for measuring end cross section. 


Figure 11. — Sheet flow of Dow-Coming 705 oil from a 3.4 cm 
wide slit. t c = 109 p.m, /Aro = 3.4, UW = 9.6, vv Q = 11 .1 m/s. 
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Figure 14. — Comparison of experimental and theoretical edge 
cross section dimensions. 
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Flattened out or cigar like cross-sectional shape 

y 



Rebuilding or peanut like cross-sectional shape 



Larger, rebuilt, nearly circular cross-sectional shape 
Figure 15. — Schematic of edge cross-section variation in flow direction. 
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(a) Edge cylinder at z = 1 .5 cm ( zJL = 0.17). 



(b) Edge cylinder at z = 2.6 cm (z/L = 0.29). 

Figure 16. — Photographs of plan edge views of sheet 
edge. W = 3.4 cm, r = 1 00 p.m, L = 8.9 cm. 
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